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摘要 : 广义 相对 论 是 继 牛 顿 万 有 引力 理论 后 ， 在 强 引力 场 、 强 动态 时 空 区 域 更 精准 的 引力 理论 ， 
大 量 的 天 文 观测 验证 了 广义 相对 论 的 正确 性 。 但 广义 相对 论 与 量子 理论 在 逻辑 体系 上 不 相 容 ， 广 
义 相对 论 还 存在 时 空 奇 点 的 理论 问题 ， 宇 宙 学 中 还 存在 暗物质 和 暗 能 量 的 引力 本 质问 题 等 。 在 使 
用 Yunes 等 人 提出 的 Einstein-aether 理论 修改 引力 波 模板 的 基础 上 ， 用 Fisher Matrix 的 参数 
古 计 法 ， 探 讨 了 几 个 典型 双星 系统 对 修改 引力 模板 的 参数 估计 所 能 达到 的 参数 估计 精度 ， 以 及 
LIGO 和 VIRGO 已 经 发 现 的 11 个 双星 并 合 事件 应 对 修改 引力 的 限制 能 力 ;， 最 后 ， 使 用 匹配 因 
子 的 方法 ， 探 讨 修改 引力 的 参数 对 引力 波 模 板 的 影响 。 如 果 Einstein-aether 理论 只 是 微小 修正 
广义 相对 论 引力 波 波形 ，11 个 双星 引力 波 事件 可 以 对 Einstein-aether 理论 的 3 个 参数 给 出 限制 ; 
如 果 Einstein-aether 理论 波形 与 广义 相对 论 引力 波 波形 存在 明显 偏离 ， 对 Einstein-aether 理论 
的 检验 需要 更 仔细 的 引力 波 波形 研究 和 数据 分 析 。 
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1 引 言 


广义 相对 论 是 继 牛 顿 万 有 引力 理论 后 ， 在 强 引 力 场 、 强 动态 时 空 区 域 更 精准 的 引力 
申论。 广义 相对 论 预言 的 很 多 物理 现象 ， 包 括 星 光 偏 折 、 引 力 红 移 、 夏 皮 罗 时 间 延 迟 、 参 
考 系 被 拖 搜 的 效应 和 黑洞 等 都 一 一 得 到 验证 。 在 广义 相对 论 发 表 100 年 后 ，2015 年 9 月 
15 日 引力 波 被 美国 的 激光 干涉 仪 引 力 波 天 文 台 (the laser interferometer gravitational-wave 
observatory, LIGO) 成 功 直接 探测 到 。 引 力 波 的 成 功 探测 完成 了 广义 相对 论 最 后 一 个 理论 预 
言 的 实验 检验 ””， 也 为 人 们 观测 天 体 和 宇宙 提供 了 一 种 全 新 手段 ””。 

然而 ， 广 义 相 对 论 仍然 存在 各 种 各 样 的 问题 。 暗 物质 和 暗 能 量 始终 是 广义 相对 论 这 片 广 
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春天 空 上 的 两 杀 乌 云 。 对 于 推动 宇宙 加 速 膨 胀 起 到 重要 作用 的 暗 能 量 ， 我 们 只 能 考虑 有 效 
场 论 ， 唯 象 层面 上 在 标准 广义 相对 论 下 加 入 一 些 修正 (通常 在 标准 广义 相对 论 下 额外 考虑 
个 标量 场 )， 最 典型 的 是 k-essence 理论 " 。 对 于 暗物质 同样 如 此 ，1933 年 ，Firtz a 
现 后 发 座 星 系 团 中 各 个 星系 团 的 速度 弥散 太 高 ， 最 开始 认为 只 是 一 些 不 发 光 的 物质 。20 世 
纪 80 年 代 左 右 ， 大 量 的 观测 表明 ， 这 种 不 发 光 的 物质 3 a E 
这 对 基础 物理 提出 了 一 个 很 大 的 挑战 ， 这 种 物质 虽然 具有 很 强 的 宏观 引力 效应 ， a 
域 上 与 物质 相互 作用 很 弱 ”。 另 外 ， 广 义 相对 论 与 量子 理论 在 逻辑 体系 上 不 相 容 ” J 
义 相对 论 自身 预言 时 空 奇 点 的 存在 ， 这 些 都 暗示 着 广义 相对 论 的 理论 缺陷 。 

对 广义 相对 论 所 有 理论 预言 的 成 功 检验 既 表明 了 广义 相对 论 的 成 功 ， 也 给 入 们 带 来 机 
会 寻找 广义 相对 论 失效 的 物理 条 件 。 牛 顿 万 有 引力 理论 适用 范围 对 应 的 弱 场 和 低速 ， 都 是 相 
对 于 广义 相对 论 黑洞 致密 程度 对 应 的 强 场 和 相对 论 预言 的 极限 光速 而 言 。 可 以 想象 ， 广 义 相 
对 论 失 效 的 物理 条 件 是 相对 于 一 种 新 理论 给 出 来 的 特征 物理 量 而 言 ， 所 以 寻找 广义 相对 论 
失效 物理 条 件 的 最 好 方法 就 是 对 比 不 同 于 广义 相对 论 的 引力 理论 ， 用 实验 数据 来 进行 理论 
选择 ”。 引 力 波 实验 正好 对 应 强 引 力 场 、 强 动态 时 空 区 域 的 物理 条 件 。 引 力 波 实验 正在 以 前 
所 未 有 的 标准 检验 着 广义 相对 论 ， 同 时 众多 的 其 他 引力 理论 被 引力 波 实验 否定 ”"””。 毫 无 
疑问 ， 科 学 研究 进展 使 人 们 对 引力 基本 物理 有 了 更 深刻 的 理解 。 

本 文 以 En 理论 为 例 ， 探 讨 引 力 波 对 超越 广义 相对 论 引 力 理 论 的 实验 检验 。 
Einstein-aether 理论 的 引力 作用 量 可 写 为 : 


4 ôo H v 
a /av (R+ MESU") ， (1) 


Mee = cl1g? “gv 二 c20567 十 036060 + elf? U” guv 3 (2) 


RH, UM 是 以 太 对 应 的 四 速 ，gu 是 时 空 度 规 ， 尺 是 里 奇 标量 ， 06 是 是 Kronecker delta 
wE, Gea 是 Einstein-acther 理论 的 有 效 引 力 常数 ， 有 别 于 牛顿 的 万 有 引力 常数 ， 
C1, C2, C3, C4 是 Einstein-aether 理论 的 4 个 参数 ， 如 果 这 4 个 参数 都 等 于 0， 则 Einstein- 
aether 理论 回 到 广义 相对 论 。Einstein-aether Hit SBA: 


Ce C= , 


C14 三 C1 十 C4 , Ci23 三 C1 十 C2 十 C3 , 


appN 8(c3 十 clc4) (3) 
2c1 — C+C- 

ppN _ OPN (e + 2e3 -cq)(2cl + Beg + c3 + ca) 

2 2 (2 B C14) C123 


脉冲 双星 的 观测 ”和 宇宙 学 的 观测 ”将 其 限制 为 : 


) 


aPPN| < 10-4 ，|agP < 10-7 , 


(4) 
全 <10? , e <103. 
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文章 余下 部 分 的 安排 如 下 : 第 2 章 介 绍 笔者 提出 的 Einstein-aether 理论 中 圆 轨 道 双 星 
系统 的 后 爱 因 斯 坦 引 力 波 波形 形式 ; 第 3 章 介绍 本 文 使 用 的 Fisher Matrix 参数 估计 方法 ; 
第 4 章 给 出 引力 波 对 Einstein-aether 理论 中 参数 的 限制 结果 和 检验 Einstein-aether 理论 的 
结果 ; 第 5 章 是 小 结 和 讨论 。 


2 Einstein-aether 理论 中 圆 轨 道 双 星系 统 的 后 爱 因 斯 坦 引 力 波 波形 


根据 文献 [中 的 计算 结果 ， 结 合 圆 轨 道 双星 系统 的 TaylorF2 引力 波 波 形 模型 ， 我 
们 把 Einstein-aether 理论 中 国 轨 道 双星 系统 的 后 爱 因 斯 坦 引力 波 波 形 写 成 如 下 形式 


h(f)=— ae Aan 5 fT tf) 4 Afi ， (5) 
1225 
T 3 1X 
= Wi(f) = 2nft. — pe — 71/44 5 [8-10 u? — Bo + — >， Ut (Dn ln lnui)| , (6) 
256u7 7 4 
= (nM f)” , (7) 
3 i 
Wo(f) = 2nft. — 2¢.— 71/44 1288 | B-n? uz? — By + 7 ten 十 1 hw) (8) 
= (nM f)? , (9) 
Po = 1 ) (10) 
lo =0 ) (11) 
Pı = 0 ’ (12) 
1, =0 ) (13) 
20 /743 11 
_ 20 (743 11 14 
S pa 2 (eon) ， (14) 
= lo =0 ’ (15) 
P3 = 一 16rr , (16) 
l3 =0 ’ (17) 
3058673 5429 617 
=1 | | IE l 
ve Gna 1008" " 144” ) Gi 
l =0 ) (19) 
38645 65 
= 1—In(xM fiso)] , 20 
ps = n (SES — Fr) (1 Inte fi) (20) 
38645 65 
b = 3n | 2> _ 22 21 
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_ 11583 231 236 531 640 2 6848y 6848 n4 


ue 4694215 680 3 21 21 
15737765635 22557? 76055 , 127825 
n’ n? (23) 
3048192 12 1728 1296 7 ° 
6848 
ham (24) 
_ (77096675 ，378515 74045. (05) 
Pr 254016 | 1512.’ 76}? 
ly = 0 ry (26) 
1 
fiso = 632nM ) (27) 
Kp, y 是 欧 拉 常数 ，M = ma + ms 是 双星 总 质量 ，n = 一 2 是 双星 系统 的 对 称 


质量 比 ，fis。 是 最 后 稳定 轨道 的 频率 ，f 是 引力 波 的 频率 ，po,p1,p2,p3, pa,ps;pe,p7 和 
lo, l1, l2, ls, la, ls, lo, l7 是 引力 波 相 位 王 M Wy 的 展开 系数 。 本 文 提出 的 波形 模型 与 文献 2 
略 有 不 同 。 文 献 四 关注 后 牛顿 近似 ， 波 形 公式 (包括 式 0, @), (2D) 和 (22) 中 的 质量 
都 是 咽 呈 质量 (chirp mass)。 我 们 关注 波形 的 准确 性 ， 使 用 了 总 质量 。 我 们 的 波形 模型 在 
a= bı = fo 二 0 时 ， 严 格 回 到 LIGO 数据 处 理 使 用 的 TaylorF2 波形 模型 。 在 唯 象 层次 上 ， 
我 们 的 上 述 模型 很 容易 直接 推广 到 LIGO 数据 分 析 软 件 库 中 最 先进 的 IMRPhenom 波形 模 
型 下 文献 [到 中 的 公式 (2) ( 即 文献 (25) 的 公式 (2.5)) 里 的 yn 和 ww 由， 以 及 文献 团 中 
WI Y- 和 wo 与 本 文采 用 的 参数 对 应 关系 如 下 : 


_ 3 2/5 -7/3 
$2 =0, Bg = Bane)" 
3 
w=", 0V = — -zbo ， 
128 
E (n—5)/3 5 (28) 
n= M) “T n zlal M ) >= , 
Ya = g MEA [on tinfed =0 
(D) 1 —(n—5)/3 (1) 
n = -z (nM In yn 一 0 . 
vl = (nm) Y 


考虑 到 GW170817 的 多 信使 天 文学 对 引力 波 速度 的 精确 测定 GERA A 10-18) 
我 们 直接 把 引力 波 速度 设 定 为 光速 。 于 是 将 上 述 波形 表达 式 与 文献 [ 的 计算 结果 相 结合 ， 
我 们 得 到 : 


a = Fsiny + YyLFr sin Y + yx1Fx cosy + yx2Fx sinw + yyıFy + 7v2Fysiny , (29) 
Bo = 0 》 (31) 


EA, y Æ Hike A, Fy, Fr, Fx, Fy 是 引力 波 探测 器 对 相应 引力 波 偏振 模式 的 响应 函 
Bl, Vy. L,Yx1,Yx2,YY1; YY2 是 Binstein-aether 理论 中 相应 引力 波 偏振 模式 的 强度 系数 。 相 


V1 


2 


06.0034 
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应 函数 和 偏振 模式 强度 系数 可 以 表达 为 Einstein-aether 理论 参数 c1, c2, c3, c4 WAM, FAK 
RIAN BES ICM [TB]. ce 和 ks 是 Einstein-aether 理论 给 出 的 参数 ， 也 是 cl cz, cas, ca 
的 函数 ， 分 别 定义 在 文献 [四 中 的 公式 (4.6) 和 (6.7)。 广 义 相 对 论 对 应 es = 0 M kg = 1。 本 
文 接 下 来 介绍 如 何 利用 引力 波 探测 来 限制 a, 8-1 和 Boo 

高 级 LIGO 探测 器 自 2015 年 9 月 开机 以 来 ， 已 经 历 了 三 次 观测 运行 ， 其 探测 灵敏 度 还 
在 不 断 提高 。 高 级 LIGO 探测 器 的 噪声 功率 谱 可 以 近似 成 下 面 的 形式 : 


So Fees 42.4 x 107622750 + 0,0817469 


1 一 0.23z2 + 0.076 42+ 
Salf) = BE ( 1 + 0.17x2 )] 


oo , f < 10 Hz 
x= f/215, So = 107" . (33) 


RP, f 为 引力 波 频 率 ，z 为 约 化 频率 ，56 为 总 的 噪声 水 平 。 


f>10Hz . (32) 


3 Fisher Matrix 的 参数 估计 方法 


引力 波 探测 器 所 测 得 的 数据 d 包含 信号 s MEE n, B d(t) = s(t) 十 n(t)。 在 理想 情况 
下 ， 探 测 器 的 噪声 是 高 斯 稳 态 噪声 ， 即 d(t) 一 s(t) 是 一 个 随机 过 程 ， 而 且 各 种 联合 概率 分 布 
都 是 与 时 间 无 关 的 高 斯 分 布 。 

基于 引力 波 探测 器 的 灵敏 度 ， 我 们 定义 两 组 时 间 序 列 h(t) 和 ho(t) 的 内 积 如 下 : 


fmax hih 


(hı (t), ho(t)) = 4R 


daf , (34) 
fmin n 

这 里 hy Al ho KER "(RAIL ER, “eR, “ROR, (fmin, fmax) 

对 应 探测 器 的 敏感 频率 范围 。 对 于 高 斯 稳 态 噪声 n(t)， 其 内 积 应 该 满足 高 斯 分 布 : 


p((n(t), n(t))) ce ne) (35) 


于 是 我 们 就 有 条 件 概率 分 布 : 


pldlh) cc ae (hah ， (36) 


对 于 引力 波 探测 ， 我 们 关心 的 是 基于 探测 数据 d(t)， 获 得 信号 h(t) 的 概率 是 多 大 ， 也 就 是 
条 件 概 率 p(h|d)。 根 据 贝 叶 斯 定理 ， 我 们 可 得 : 
dp 
UA gap 
上 式 分 母 实 为 分 子 对 应 的 概率 归 一 化 因子 。 在 贝 叶 斯 统计 学 里 ，p(h|q) 被 称 为 后 验 概率 ， 
p(h) 被 称 为 先 验 概率 ，p(dlh) 被 称 为 似 然 函数 。 


(37) 
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对 于 某 一 种 类 引力 波源 ， 如 双星 系统 ， 引 力 波 波形 除了 是 时 间 的 函数 外 ， 还 依赖 于 波源 
的 参数 ， 如 双 黑 洞 的 质量 、 自 旋 和 轨道 离心 率 等 。 我 们 把 这 些 参数 记 作 和 ， 则 可 得 AC, t) 
概率 分 布 p(h|d) 和 p(h) 是 相对 于 入 的 函数 ， 即 是 参数 空间 上 的 函数 。[ p(djh)p(h) 是 对 整 
个 参数 空间 积分 。 

后 验 概率 plhd) 作为 参数 空间 上 的 函数 ， 其 最 大 值 所 对 应 的 点 被 称 为 最 可 几 参 数 和 0。 
eee No 对 应 最 好 匹配 的 那 组 参数 二 。 我 们 再 进一步 假设 

号 很 强 ， 即 信 品 比 很 高 ， 品 声 n 基本 可 忽略 不 计 ， 则 有 qa = A(X), (n,n) & (n,h) & 
J ) a (n, OO wh) 守 0。 如 果 我 们 简化 认为 先 验 均匀 分 布 p(h) x 1， 则 根据 式 (Bo) 和 
(87) 可 得 : 


ph) = Ce ee? ， (38) 


其 中 ，C 是 概率 分 布 的 归 一 化 因子 。 我 们 考虑 Inp(h(X)|q) 在 Xo 点 处 的 泰勒 展开 近似 : 
In p(h(A)|d) = nC — (d — h,d — h) /2 ~ In C + A — NOih, ðh) (A — NM) . (39) 


上 式 右边 的 矩阵 (Och, Ash) 刚好 对 应 统计 学 里 的 Fisher Matrix。 也 就 是 说 在 高 信 噪 比 
条 件 下 ， 后 验 概率 分 布 在 最 可 几 参 数 点 可 近似 为 协 方差 矩阵 X F Fisher Matrix 逆 
(Oxih, Osh) = Ez HW. 

本 文 接 下 来 使 用 Fisher Matrix 的 方法 来 估计 引力 波 测量 的 参数 误差 


AX = V5 . (40) 


本 文 使 用 的 Einstein-aether 理论 中 ， 圆 轨道 双星 系统 的 后 爱 因 斯 坦 引力 波 波 形 模型 包 
含 参 数 (A, M,N, te, Qc, a, p_i6o)， 它 们 分 别 是 引力 波幅 度 (与 引力 波源 的 光度 距离 和 方位 
角 等 有 关 )、 咽 嗽 质量 、 对 称 质 量 比 、 并 合 时 间 与 相位 ， 以 及 Einstein-aether 理论 对 应 的 后 
爱 因 斯 坦 参 数 。 我 们 在 计算 中 发 现 A 与 其 他 参数 的 关联 度 很 弱 ， 该 参数 是 否 参 与 Fisher 
Matrix， 所 得 其 他 参数 的 误差 基本 相同 。 使 用 高 级 LIGO 探测 器 的 灵敏 度 ， 对 典型 双星 系统 
参数 误差 结果 如 表 四 所 示 。 我 们 的 结果 与 文献 2 的 一 致 。 由 于 我 们 的 波形 模型 与 其 的 细微 
差别 导致 了 定量 上 的 细微 差别 ， 从 差别 可 以 看 出 ， 我 们 采用 更 精确 的 波形 模型 (等 效 到 3.5 
阶 后 牛顿 ) 求 出 的 参数 误差 比 采用 文献 2 的 波形 模型 (等 效 到 2 阶 后 牛顿 ) 算得 的 略 小 。 


表 1 使 用 Fisher Matrix 方法 针对 高 级 LIGO 灵敏 度 对 典型 双星 系统 进行 参数 估计 的 结果 
mi m2 Ate Ade AM/M An/n Aa Ab-ı Abo 
/ms 
14 1.4 0.331 0.925 3.84x107? 1.64x107? 0.547 2.79x107ë 0.262 
2.0 1.0 0.335 0.976 3.86x 107? 1.40x107? 0.560 3.47x107 0.296 
1.8 1.2 0.332 0.936 3.98x107? 1.52x107? 0.552 3.29x10-5 0.283 
1.4 1.5 0.331 0.921 3.94x107? 1.61x10-2 0.547 2.99x10-5 0.269 


注 : 假设 信 噪 比 等 于 10， 质 量 单位 是 Moe 
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KF te 和 内 这 两 个 参数 ， 我 们 还 尝试 了 使 用 下 统计 方法 来 处 理 ”。 对 于 参数 误差 ， 使 
JF 统计 与 否 ， Binstein-aether 理论 对 应 的 后 爱 因 斯 坦 参 数 限制 结果 基本 没有 差别 。 在 下 
文 计算 中 ， 我 们 得 到 的 结果 都 是 没有 使 用 下 统计 的 计算 结果 。 


4 对 Einstein-aether 理论 进行 限制 


高 级 LIGO 主导 的 LIGO 科学 合作 组 织 第 一 次 观测 运行 从 2015 Æ 9 H 12 日 到 2016 年 
1 月 19 日 ， 共 发 现 GW150914, GW151012 和 GW151226 三 个 双 黑 洞 并 合 事件 。 高 级 LIGO 
主导 的 LIGO 科学 合作 组 织 第 二 次 观测 运行 从 2016 年 11 月 30 日 开始 到 2017 年 8 月 25 
日 结束 。 其 中 高 级 VIRGO 从 2017 年 8 月 1 日 开始 加 入 第 二 次 观测 运行 。LIGO 科学 合作 
组 织 第 二 次 观测 运行 共 发 现 GW170104, GW170608, GW170729, GW170809, GW170814, 
GW170817, GW170818 和 GW170823 八 个 双星 并 合 事件 。 

我 们 把 LIGO 科学 合作 组 织 第 一 次 和 第 二 次 观测 运行 发 现 的 11 个 双星 并 合 事件 对 
Einstein-aether 理论 进行 参数 估计 的 结果 列 在 了 表 回 中 。 因 为 只 有 GstLAL 数据 处 理 流水 线 
发 现 了 所 有 的 11 个 事件 ， 所 以 表 回 里 所 列 信 噪 比 SNR 和 双星 系统 的 质量 m 和 me 都 是 由 
GstLAL 数据 处 理 流水 线 所 得 ~”。 我 们 将 LIGO-VIRGO 第 一 次 和 第 二 次 探测 到 的 11 个 引 
力 波 事件 的 参数 代入 到 式 回 一 四 中 ， 得 出 带 有 Einstein-aether 理论 的 3 个 参数 的 模拟 波 
6, HEH Fisher Matrix 方法 估计 引力 波 事件 对 这 3 个 参数 的 限制 能 力 。 从 该 表 结 果 可 以 
看 出 ， 双 星系 统 的 质量 越 小 ， 则 对 Einstein-aether 理论 的 参数 限制 越 严格 ;这 是 因为 质量 
越 小 ， 落 到 高 级 LIGO 和 高 级 VIRGO 探测 器 频段 中 的 引力 波 信号 越 多 ， 对 Einstein-aether 
里 论 的 验证 就 越 有 效 。 当 双星 系统 的 质量 增 大 ， 直 到 并 合 频率 落 到 探测 器 频段 中 ， 虽 然 并 合 
信和 号 会 较 强 ， 但 我 们 所 使 用 的 引力 波 波形 模型 无 法 利用 上 并 合 短信 号 ， 导 致 所 得 限制 不 强 。 


= 


#2 11 个 双星 并 合 事件 对 Einstein-aether 理论 进行 参数 估计 的 结果 
ID SNR m mea Aa AB-1 ABo 
/Mo /Mo 
GW150914 24.4 35.6 30.6 0.191 1.95x10-? 0.392 
GW151012 10.0 23.2 13.6 0.510 7.77x10-3 0.489 
GW151226 13.1 13.7 7.7 0.407 1.26 x 107 0.351 
GW170104 13.0 30.8 20.0 0.376 1.59x10-? 0.467 
GW170608 14.9 11.0 7.6 0.358 7.48x1074 0.293 
GW170729 10.8 50.2 34.0 0.414 9.96 x 107? 1.64 
GW170809 12.4 35.0 23.8 0.385 2.65107? 0.609 
GW170814 15.9 306 25.2 0.301 1.73x10-? 0.430 
GW170817 33.0 1.46 1.27 0.166 8.07x107 0.0779 
GW170818 11.3 35.4 26.7 0.418 3.44107? 0.739 
GW170823 11.5 39.5 29.0 0.404 467x107? 0.901 
TE: 后 面 3 列 参数 的 估计 结果 由 我 们 的 Fisher Matrix 分 析 所 得 。 
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综合 LIGO 科学 合作 组 织 第 一 次 和 第 二 次 观测 运行 发 现 的 11 个 双星 并 合 事 件 对 
Einstein-aether 理论 进行 参数 限制 的 结果 ， 可 以 得 到 : 


a<0.166, Bı <807x10°°, 6) <7.79x10- 2 . 


但 我 们 要 注意 ， 以 上 限制 结果 基于 假定 广义 相对 论 的 引力 波 模板 可 以 完全 匹配 探测 器 所 测 
得 的 信号 。 引 力 波 对 a 参数 限制 不 好 ， 对 8_1 和 Bo 限制 不 错 ， 这 是 因为 a 与 引力 波 的 幅度 
相关 ，B6_1 和 By 与 引力 波 的 相位 相关 。 而 匹配 滤波 本 来 就 只 是 敏感 依赖 于 相位 。 

为 了 估计 广义 相对 论 波形 模板 匹配 Einstein-aether 理论 引力 波 信号 的 程度 ， 我 们 引入 
引力 波 匹配 因子 FETT., AEA RR h(t) 和 ho(t)， 它 们 间 的 匹配 因子 为 ; 


(hi(t), ho(t)) 
F | y 
(h(t), hit) - (halt), halt) (41) 


当 用 一 个 波形 近似 另 一 个 波形 信号 时 ， 信 噪 比 会 被 降低 PP. WREERE, KEERA 
近似 就 会 失效 。 

我 们 以 GW170817 的 双星 参数 为 例 ， 研 究 广义 相对 论 与 Einstein-aether 理论 的 引力 波 
波形 匹配 情况 。 鉴 于 匹配 滤波 与 引力 波幅 度 不 是 敏感 依赖 ， 我 们 只 讨论 匹配 因子 与 6_1 
和 bo WKA, BROAD 从 图 加 我 们 可 以 看 出 ， 只 要 bM Bo 偏离 0 分别 在 10-6 和 
10-3 的 量 级 ，Einstein-aether 理论 的 引力 波 波形 与 广义 相对 论 的 波形 完全 不 匹配 。 所 以 只 
要 GW170817 的 引力 波 信号 是 广义 相对 论 描述 的 波形 ，Einstein-aether 理论 的 86_1 和 Bo 参 
数 就 可 以 得 到 很 好 的 限 和 


SSS 


=. 
o 


XxX 10° 


a) b) 


1 ”类 似 GW170817 的 广义 相对 论 引力 波 波形 和 Einstein-aether 理论 的 引力 波 波形 匹配 情况 


这 里 还 存在 一 个 问题 : 是否 GW170817 实 为 Einstein-aether 理论 描述 ， 只 是 因为 b1 
和 By 与 咽 嗽 质量 的 参数 简 并 使 得 广义 相对 论 模 板 没 有 遗漏 这 个 信号 ? 接 下 来 我 们 研究 这 个 
问题 。 为 此 ， 我 们 设 定 6_1 = 6 x 107, By = 0 M 81 = 0, Bp = 0.002， 然 后 改变 咽 呈 质 量 ， 
研究 波形 匹配 因子 ， 结 果 见 图 色 从 图 加 可 以 看 出 ， 的 确 可 以 通过 改变 咽 嗽 质量 来 使 得 广义 
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相对 论 的 引力 波 波形 与 Binstein-aether 理论 的 引力 波 波形 相 匹配 ， 即 8_1 和 By 与 咽 哆 质量 
存在 参数 简 并 。 所 以 在 不 确定 所 探测 到 的 引力 波 信号 可 由 广义 相对 论 描述 的 情况 下 ， 引 力 波 
的 探测 不 能 被 认 已 经 对 Einstein-aether 理论 给 出 了 很 好 的 限制 ， 现 在 我 们 只 能 说 广义 相对 
论 能 很 好 地 解释 所 有 探测 到 的 引力 波 信号 ”。 


1.0 
0.8 
= 06 
0.4 
0.2 : : : 1 : 0.2 
2.796 2.798 2.800 2.802 2.804 2.796 2.798 2.800 2.802 2.804 

M, M, 
a) b) 

YE: 这 里 我 们 使 用 了 n = 0.25 和 A = 2.8Mo 作为 参考 。a) 比较 了 6_1 = 6 x 1077, 8o = OF 


Einstein-aether 理论 与 不 同 咽 嗽 质量 的 广义 相对 论 引力 波 波形 匹配 情况 ; b) 比较 了 6_1 = 0, Bo = 0.002 时 
Einstein-aether 理论 与 不 同 咽 吹 质量 的 广义 相对 论 引力 波 波形 匹配 情况 。 


2 “Einstein-aether 理论 中 6_1 和 Bo 与 咽 嗽 质量 的 参数 简 并 情况 


5 结论 和 讨论 
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自 2015 年 LIGO 首次 直接 探测 到 引力 波 以 来 ， 


三 次 观测 ， 共 
引力 波 天 文学 可 以 用 来 研究 基本 物理 学 中 的 引力 相互 作用 一 一。 


功 应 用 到 引力 理论 的 检验 中 ， 若 干 不 同 于 广义 相对 论 的 引力 理论 得 到 了 很 强 的 P 


高 级 LIGO 和 高 级 VIRGO 一 起 进行 了 


探测 到 几 十 例 引力 波 事件 ， 包 括 双 黑洞 、 双 中 子 星 和 黑洞 -中 子 双星 系统 的 并 
合 ， 以 引力 波 为 观测 手段 的 引力 波 天 文学 被 很 好 地 建立 起 来 。 


在 文献 R 的 计算 结果 基础 上 ， 结 合 TaylorF2 引力 波 波 形 的 模式 ， 


i=! 


GW170817 就 被 成 


民 制 |。 
我 们 提出 了 


Einstein-aether 理论 中 双星 系统 的 引力 波 波 形 模型 。 基 于 所 提出 的 波形 模型 ， 我 们 利用 
LIGO 科学 合作 组 织 公 开 的 第 一 、 第 二 次 观测 运行 结果 数据 中 的 


Fisher Matrix 方法 ， 结 合 


11 个 双星 系统 ， 对 Einstein-aether 理论 进行 了 限制 、 研 究 。Fisher Matrix 方法 只 对 参数 微 
小 偏离 ， 或 者 说 只 在 Einstein-aether 理论 微小 修正 广义 相对 论 引 力 波 波形 的 | 


这 种 情形 下 ， 使 用 
与 广义 相对 论 相似 的 结果 ， 所 以 我 们 只 需要 使 月 


Einstei 


n-aether 理论 的 引力 波 波形 模板 对 引力 波 实际 数据 的 处 理会 得 出 
H Fisher Matrix 方法 对 后 验 概率 分 布 做 微 


青 形 下 有 效 。 


扰 研 究 即 可 。 这 样 的 条 件 使 得 我 们 可 以 把 LIGO 已 经 测 出 的 11 个 双星 对 应 的 参数 输入 到 


Fisher Matrix 方法 进行 分 机。 反之 ， 


且 参 数 微小 偏离 不 成 立 ，Fisher Matrix 


方法 将 失效 。 
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我 们 采用 匹配 滤波 中 匹配 因子 的 方法 ， 探 讨 广 义 相 对 论 引 力 波 波形 与 Einstein-aether 理论 
的 引力 波 波形 存在 显著 偏离 的 可 能 性 及 其 对 引力 波 数据 处 理 的 影响 。 

结果 表明 ， 如 果 Einstein-aether 理论 只 是 广义 相对 论 的 微小 偏离 (在 引力 波 波形 意 
义 上 )， 则 我 们 可 以 把 Einstein-aether 理论 中 的 参数 限制 为 a < 0.166, 8 < 8.07 x 
1076, Bo < 7.79 x 10-?。 但 只 是 基于 引力 波 探测 ， 我 们 还 不 能 完全 排除 Einstein-aether # 
论 显著 偏离 广义 相对 论 的 可 能 性 ， 这 时 Einstein-aether 理论 的 参数 和 双星 系统 的 咽 喇 质量 
存在 简 并 。 我 们 计划 在 后 续 工作 中 ， 使 用 贝 叶 斯 模型 选择 的 方法 来 研究 这 种 简 并 的 可 能 性 有 
多 大 。 
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On the Test Ability of Einstein-aether Theory with 
Gravitational Wave Events in LIGO 01/02 


LIU Gang-qiang, CAO Zhou-jian 


(The Department of Astronomy, Beijing Normal University, Beijing 100875, China) 


Abstract: General relativity is a more accurate gravitational theory for strong gravitational 
field and dynamic space-time region after Newton’s theory of gravity. A large number of 
astronomical observations have verified the theory of general relativity. However, general 
relativity and quantum theory are not compatible logically. There are also some theoretical 
problems in general relativity, such as the singularity of space-time, the problem of dark 
matter and dark energy in cosmology. Based on modified waveform of the Einstein-ather 
theory proposed by Yunes and his coworkers, Fisher matrix method is used to discuss the 
parameter estimation accuracy of several typical binary systems, and the ability of 11 binary 
merging events discovered by LIGO and Virgo are used to constrain the modified gravity. 
Finally, using the matching factor to discuss the influence of the modified gravitational 
parameters on the gravitational wave template. If Einstein-ather theory just slightly modifies 
the gravitational wave waveform of general relativity, 11 binary gravitational wave events 
can limit the three parameters of Einstein-ather theory. If the waveforms of Einstein-ather 
theory deviate significantly from the waveforms of general relativity, the test of Einstein- 


ather theory needs more careful study of gravitational wave waveforms and data analysis. 


Key words: gravitational wave; general relativity; Einstein-ether theory 


